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1.  Introduction 

Several  test  statistics  in  multivariate  analysis  are 
based  upon  certain  functions  of  the  elements  of  the  sample 
covariance  matrix  or  sample  correlation  matrix.  For  example, 
the  sample  correlation  coefficient,  partial  and  multiple 
correlation  coefficients,  and  various  transformations  of 
the  sample  correlation  coefficients  depend  upon  the  elements 
of  the  sample  covariance  matrix  and  correlation  matrix.  The 
exact  distributions  of  many  of  these  statistics  are  quite 
complicated  and  so  there  is  a  need  to  derive  asymptotic 
expressions  of  the  above  functions.  In  this  paper,  we 
consider  asymptotic  joint  distributions  of  functions  of  the 
elements  of  the  noncentral  Wishart  matrix  and  the  associated 
noncentral  correlation  matrix. 

In  Sectioii  2  of  this  paper,  we  derive  the  asymptotic 
joint  distribution  of  certain  functions  of  the  elements  of 
the  noncentral  Wishart  matrix.  The  first  term  in  the  asymp¬ 
totic  expression  is  the  multivariate  normal  density,  whereas 
the  second  term  involves  partial  derivatives  of  the  multi¬ 
variate  normal  density.  Similar  expressions  are  given  for 
the  case  of  the  noncentral  correlation  matrix.  The  method 
used  involves  expanding  the  functions  in  terms  of  Taylor  series 
and  computing  the  characteristic  functions.  Olkin  and 
Siotatii  (1976)  obtained  the  first  term  in  the  asymptotic 
joint  distribution  of  functions  of  the  elements  of  the  central 
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correlation  matrix.  Siotani  and  Hayakawa  (1964)  obtained 
the  first  terms  in  the  asymptotic  joint  distributions  of 
the  partial  and  multiple  correlation  coefficients  by  express¬ 
ing  them  as  functions  of  the  elements  of  the  central  Wishart 
matrix.  Konishi  (1979)  obtained  the  first  two  terms  in  the 


asymptotic  joint  distribution  of  various  functions  of  the 
central  correlation  matrix;  these  results  are  special  cases 
of  the  results  given  in  this  paper.  In  Section  3  of  this 
paper,  we  studied  the  accuracy  of  the  asymptotic  expressions 
given  in  Section  2  for  some  special  cases.  Asymptotic  ex¬ 
pressions  for  the  joint  distributions  of  functions  of  the 
elements  of  the  sample  covariance  matrix  are  given  in  Section 
4  when  the  underlying  distribution  is  not  multivariate  normal. 
Finally,  applications  of  the  results  of  this  paper  are  dis¬ 


cussed  in  Section  5. 
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2.  Joint  Distribution  of  the  Functions  of  the  Elements 


of  the  Noncentral  Wishart  Matrix 


Let  X-j . be  distributed  independently  as  multi¬ 
variate  normal  with  mean  vectors  pc . and  covariance 

-1  n 

matrix  £  =  (a..).  Then,  the  distribution  of  S =  7  X.X'=(S..) 

ij  ji'j-J-j 

is  known  to  be  the  noncentral  Wishart  matrix  with  n  degrees 


of  freedom  and  E(S/n)  =  E  +  (M/n)  =  Q,  with  noncentrality 
n 

parameter  M=  ][  y  . p ’  =  n(  v .  .  )  .  Now,  let 


j  =  l 


“J~J 


ij 


T.(S/n)=  T.(s,,,...,s  ,  s10I . , . ,s,  ,s„. . s  n  ) 

1  '  l  11’  ’  pp’  12’  ’  lp’  23  ’  p-l,p 

for  i  =  l,...,k,  where  s..  =  S../n.are  analytic  functions  in 

...  i.r 

the  neighborhood  of  (1  =  (ok  .).  Also,  let 


J1J2 


1+6 


jlJ2 


9  s 


J1J2 


T  •  ( S/n ) 


(S/n)  =  n 


(2.1) 


t(i) 

j 1^2  ^3^4 


1  +  6  . 


J3J4 


1  +  6  . 


Jlll2 


9s  .  .  9s  . 


.i3j4 


J1J2 


T.(S/n) 


(S/n)=ft 


where  6,  ,  is  given  by 
hk 


’hk 


1  h=k 


n  h^k. 


The  Taylor  expansion  of  T^(S/n)  about  SI  is 


Ti(S/n) 


P  P 
T.(Q)  +  l  l 

h=1  V1 


.(i) 


jjj2  ((Sj]J2/n)“<i)JiJ2) 


♦i 


P 

.1  .  .  « 

J 1  *  ^  2  ’  ^  3  ’  ^  4 


,(i) 


WW<SW''>^lJ2><<W*)'*V< 


+  higher  order  terms 


(2.2) 


P 

where  }  denotes  the  summation  over  all 

jj. j2>  * • • > Ju 

values  of  jj , jg* • • • * Ju  varying  from  1  to  p. 

Now,  let 

I,.  =  /in  (TjCS/n)  -  T.(il)  }.  (2.3) 


3  A 


w 


Using  (2.2),  we  obtain  the  following  expression  for  the 

joint  characteristic  function  of  L1  =  (L^ . L^): 

k 

i/>(t)  =  E{exp(  i  T  t  .L.)> 

i=l  0  J  (2.4) 

=  Ej ( t )  +  E2(t)  +  0( n-1 ) 
where  t'  =  (t^,...,t^),  and 


B=  l  t.(u('j  )  =  l  lj  A 

i=l  1  J1J2  i= I 


(D 


K  (  t )  *  E 

2  - 


^  S  ?  VilLi  i((Si  i  /n)~wi  1  )((si  i  /n)-wi 
^  t-I  J1.J2»j3*J4i  Jl  * 12  l^2  Jgj4  J 


(2.6) 


exp 


[i/S  l  1  t.a' Y  («  )]J 

il  =  1  jj,j2  12  12  1  2  i 


Starting  from  (2.4),  we  obtain  the  following  asymptotic 
expression  for  the  joint  characteristic  function  of  ,... : 
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nt)  «  «  xr>(  i ' (j  t.  >  n  +-1—  y  i  t .  h, 

:Vn  i  =  i  1 


(2.7) 


,  k 

+  —  }  i  L  t.  t.  (h9  +  h„)  +  (n 

/ifll*i2*i3  1  2  3 


where 


(U)  09)  (i.)  (i.,) 

Q  =  (Qj  :  )  .  Q,  j  2  tr  It  R  “  Htr  H  i|»  ~ 

J  |  1 |  1 2 


I*  I*  ).  a;  i  .  i  i  (  ,J  i  i  ,0  i  !  f(>i  i  a)  i  j  +0  .  .  V.  .  t  ill  . 

.1  ,  .  J2  '  i-t  ’  i  1  2  '  ^  1:5  '  2'  1  '  2‘  2  •*  2,l3  -1  |  ■'.!  •12J-1  J  I 


a  (  i  ,  )  (  i9)  (  i.. )  (  i ,  )  (  i9)  (  i .,) 

h.jj-  ~  lr  |{  It  ~  R  +■  4  tr  It  R  ij> 


(2.H) 


l»  (  i  ,  )  (  i  9 )  (  i ,, )  (  i )  ( i ., ) 

h  -  2  )  a  :  /  \*  +  T  f  )  ( U  .  +  T  ) 

3  .1,  .  j2,.i.r.i1  ,ilJ2*J:tJ4  Jl^2  J3J4  j.vM 


R(i)  =  A<i)E  >  =  A(l)  M/n, 


=  E  A^ft  ,  1  ^  ^  A^x)E 


and  j  denotes  the  (i,j)th  element  of  matrix  U  =  (U^). 

By  inverting  the  above  characteristic  function,  we  obtain 
the  following  expression  for  the  joint  density  of  Lj,...,Lr: 


(2.9) 


.  K 

f  (  L,  ,  .  .  .  ,  L.  )  =  N(L,Q)  U  +  -7:  l  H  (L)  h 
1  K  2/n  i=l 


i1,i2,i3  123 


where 


N(L,Q)  = 


(/5?)  '  |Q|a 


(-  |  ij'Q-1^) 


(2.10) 


^..•VLWL.Q)  =  (-1)U  ^ --3-  NCL.Q) 


_  j, 

The  correlation  matrix  is  given  by  U  =  S  “S  S  11  =  (r.  .) 

0  0  ij 


where  S  =  diag . ( S  .  .  .  , S  )  .  Now,  let 


°j<R)  =  Gj(r12’r13 . rlp,r23’  •  *  *  ,r2p"  •  •  ’rp-l,p) 


(2.11) 


be  an  analytic  function  in  the  neighborhood  of 

^  fi  ^  =  P*  =  (P*j)  where  ^  *  diag.((»>jj . ^pp^'  **  we 

denote  r.  .  by  T.  .  (S/n),  Eq.  (2.11)  can  be  written  as 

K1K2  K1K2 


G j ( R )  =  Gj(T12(S/n), . . . , Tp_1  p( S/n) ) =  (G  0  T)j(S/n).  (2.12) 


Now,  let 


1 


(j)  , 

°k  k  =  SF - (;j(k°  I 

12  kk  i  K  -  P* 


Jj) 

k, k  »k.,k,  llr.  .  Sr. 


i  2  ":r  \  k,k.  k .  k. 


R  "  •>* 


in  Kq.  (2. 1  )■ 

(.)) 


1 J2 


Then  we  obtain 
I  +  -5 

V  2 


a 


,  ,  o  K  .  . 

\  J  ,  •)  2 


\(G«T)  ^(S/n) 


(R/n)=f2  (2.13) 


‘  *  ’ 1  ■  •' •< 


<j>  t(k,k2). 

1  ' 


(  j) 


(1+6,  ,  )  (1+6,  ,  ) 


J3J4 


•11J2 


'  Ji-'- 


(Gi»T)  j(S/n) 


(S/n)=a 


(J) 


( kl k2 )  r^k3k4) 


V  v  c»  c  ^  c  ° 

*  <k  ,  k  i k 2 • k 3 k 4  ’  *i  I  * 2  V 

(.1)  •  (k,k2) 


kl"k2  k3  ”4 


r,  .  . 


k.<k0  klk2.  *  1  '*  2  '  ^  3'*-1 

i  (j 


(2.14) 


where 


2  /*.  ,  »*),  , 

^  j  ^  I  ^  *  j  ^  •) 


•*1“  kl  1  J2  “  k2  °r  ~  k2’  j2  "  kl 


(klk2}  J 

J|.*3  =  > 

~  *’k  ,  k. 

•'  i  -  h  ki 


(2.15) 


•j  !  -  -'2  =  k2 


V  0 


o t  herwi  s<* 


7 


I 


Uik2) 


r 


l 

1  "’it  k 

II 

3/ 2 
k2k2 

- 1 

1  3/2 

4  ul.  , 
kik] 

U>k2k2 

<’k  k 

k  l  2 

Jl  Hi,  , 
kThl 

“>k2k2 

*i,,R  ,K2 

•1u>, 


k  .  k  . 
1  1 


for  any  3  values  of  equal 


for  Wk2 

or  Jl=i2=k2'-i3=j4  =  kl 


for  all  j  l=.i2=r-j3=-j4=kj  ’  i  = 


=  1 .2 


(j 


otherwise 


Substituting  (2.13)  and  (2.14)  in  Eq.  (2.9)  we  obtain  the 
nsyrnpl  ol  i  r  density  for  functions  of  the  elements  of  corre¬ 
lation  matrix. 


I  o  kv 


to  kj 
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3.  An  Empirical  Study  on  the  Accuracy  of  tho  Asymptotic 

Expressions 


In  this  section,  wo  study  the  accuracy  of  the 
asymptotic  expressions  derived  in  Section  !„  for 
some  special  cases. 

We  will  first  study  the  accuracy  of  the  approximation 
for  the  distribution  of  the  noncentral  chi-square.  Let 
yi  =  Sii^°ij  r°r  i  =  J .2 - -  - •  >P*  Then  y^  is  distributed  as 
the  noncentral  chi-square  distribution  with  n  degrees  of 
freedom  and  with  the  nonr.cn tral  1  ty  parameter  y  =  nv.^/o^j. 
Now  let 


exp( - y/2 )  }  — 

j=0  .1! 


2(n/^{+:>.j1.((n/2)+j) 


v(n/2)+j-l 

o 


exp(-x/2)dx  = 
(3.1) 


and  u  is  given  by 


1 

2n^2r(n/2) 


fU 

exp(-x/2) 

Jo 


(n/2)-l 


dx 


0 . 95 . 


(3.2) 


The  left-side  of  (3.1)  is  tho  probability  integral  of  yj . 

Also  the  left-side  of  (3.1)  is  equivalent  to  the  left-side  of 
(3.2)  when  y=0.  Table  1  given  below  compares  the  exact 
values  of  3  with  the  corresponding  values  obtained  by  using 
the  asymptotic  expression  (2.9). 

TABLE  1 

Comparison  of  the  Asymptotic  Expression  with  Exact 
Expression  for  the  Noncentral  Chi-square  Distribution 
n _ u _ j _ 0(  1  ) _ 0(n  °)  0(l)40(n  "^Exac.t 

5  11.071  0.  .972(5  -.0338  .9308  .95 

16.47  .1103  -.0091  .1072  .10 
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TABI.K  1  (Continued) 


n 

li 

Y 

0(1  ) 

0(n  -)  0(])+0(n“*) 

Exact 

2.67 

.  7727 

.  02-1 5) 

.75)72 

o 

oc 

10 

18 . 307 

0. 

.  5)68  1 

0260 

.1)424 

.  95 

20.53 

.1132 

- . 0082 

.  1050 

.10 

3.71 

.  7820 

.  0155) 

.7979 

.80 

20 

21.41 

0 

.  5)6-1 1 

-.0186 

.  94  58 

.95 

26.13 

.  1  1  0-1 

-.0071 

.  1  033 

.10 

5.18 

.7880 

.0105 

.7985 

00 

o 

lit) 

19 .779 

0. 

.  5)623 

-  .  015.3 

.  9470 

.  95 

30 . 38 

.  1  081) 

- . 0064 

.  1  025 

.30 

6 . 3 1 

.  75)06 

.  0083 

.7988 

.80 

The 

values 

in  the  column  ”0(1)"  give  the  values  of 

8  when 

t  he 

f  i  rs  i 

t  erm  1 n  the 

asymptotic  expression  (2.9)  is  taken 

wlier 

cas  the  column  "0(l)  +  0(n* 

)"  gives  the 

values  of 

6  when 

the 

first 

two  terms  in 

(2.9) 

are  taken. 

"ho  column 

"0( n~ " 

Hives  tiie 

cont  rl but  ion 

of  the 

second  form 

in  (2.9) 

to  the 

value  of  H 

.  The  exact 

val  vies 

gi  von  i n  <  hr 

•  last'  column  are 

Inl'cn  from  Owen  (15)62). 

'Vo  will  compare  tho  asymptotic  expansion  given  by  (2.9) 
for  tho  distribution  of  y  ^  +  y with  the  cor  respond  i  ng  exact 


expression.  When  p 


12 


0  the  distribution  of  yj+y^  1  lie 


noneentral  chi -square  with  2n  degrees  of  freedom  and  with 
Y  ]  2  !IS  ^*1<“  noncentrality  parameter  where  ■,  ^  =  nv1j/n-|j  and 

>2  =  n  v22  /°22 '  w,inn  P]  2  f  °>  t,lfl  distribution  of  y^+yg  i  y 
the  same  as  tho  distribution  of  a  quadratic  form  in  the 
nonc.entral  case.  Now,  1  ot. 


PI 


y,  ■*  y2  1 I  p.  y-|  -  y?.  1  =  p 


(3.9) 


10 


where 


Ply,  +  y2  I  U !p-0,y1=y2=0  I  =  0.95  (3,4) 

In  Table  2,  the  entries  in  the  column  "0(  1  )+0( n~~ ) "  represent 
the  approximate  value  of  6  when  the  first  two  terms  in  the 
asymptotic  expression  (2.9)  are  taken.  The  entries  in 
the  columns  "0(J)"  and  "0(n  “)*’  represent  respectively  the 
contribution  of  the  first  term  and  the  second  term  of  the 
above  approximate  value  of  8.  The  entries  in  the  column  "exact 
values"  are  computed  by  Monte  Carlo  methods  using  the  1MSL  sub¬ 
routine  GGNSM.  In  computing  the?  simulated  values,  50C0  trials 
are  performed. 


TAUI..K  2 

Comparison  of  the  Asymptotic  Expression  with  Exact  Expression 
for  the  Distribution  of  Sum  of  Correlat'd  Chi-Square  Variables 


nun 

y2 

()(  1  ) 

0  (  n  "  ) 

0(1  )+o(n*) 

Exact 

10  31.41  0. 

0. 

0. 

.  9644 

-.0186 

.  9458 

.95 

.5 

0. 

0. 

.9467 

-.0230 

.9237 

,9364 

.5 

4. 

2. 

.  7353 

.0238 

.7591 

,7574 

.9 

0. 

o. 

.9100 

-.0192 

.8908 

,9032 

.9 

4  . 

2. 

.7092 

.0310 

.7  132 

,7362 

25  67.505  0. 

0. 

0. 

.9600 

-.01 19 

.  948  1 

.95 

.2 

0. 

0. 

.9570 

-.0125 

.0444 

,9454 

.8 

0. 

0. 

.9142 

-.0126 

.9016 

,9018 

.5 

4, 

2. 

.  5730 

.  0289 

.  6019 

,  5960 

We  now  discuss 

the 

distribution  of 

the  ratio  Fq  = 

yl/y2' 

The  distribution 

of 

F0 

is  k n own 

(Bose  (1935)  for  y,=Y2: 

=0  to  be 

2n(l-p2)n/2n  J(n  +  1 )  |  f"-’(J+F2  ) 

°'  /?  F (  } n  )  I  l+d  +F2)2  -  4  (i2Fq  -|(n+1)/^ 


*(Fn) = 


(3.5) 
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Finney  (1938)  .showed  lli.it 


Pt*lF(>  i  I2)  -  1  -  I  (Sn.in)  C».U) 

where 

x  =  *  N  -  - i  r, - 75— t  1 

ut+i.  1y-  4p“j:: 

Tx(a-h)  =  Bnrrnr  f0  y,a_1  (1-y)b_1  dy- 

Krishnaiah  el.  al  .  (1900)  obtained  an  alternative  expression 
for  l*r  (  Fq  0  I  when  y ^  and  also  gave  tables  for  this 

distribution.  fn  Table  3,  we  give  the  values  of 
1*1*1  Fo  £  u  |  p ,  Y  |  ,  Y.)  I  -  H  by  using  the  asymptotic  expressions 

(2.9),  where  u  is  the  TtV.?.  critical  value  of  the  central  F 
distribution  with  (ti.nl  decrees  <.f  freedom-  Tn  the  lu.sl 
column,  the  entries  with  *  are  obtained  by  using  the  formula 
(9.0)  and  the  remaining  entries  in  this  column  are  obtained 
by  si  mill  at  ion  for  0000  trials. 

TABLE  3 

Comparison  of  the  Asymptotic  Expression  with  Exact 
Expression  for  the  Distribution  of  the  Wat.io  of 
Correlated  Chi-square  Variables 


n  u 

i> 

*1 

y2 

0(  1  ) 

0(n  *') 

0( 1 )+0( n~“  ) 

Exact 

24  1 .32  14 

0. 

0. 

0. 

.7844 

-  .0370 

.7474 

* 

.75 

.2 

0. 

0  . 

.  7892 

.0373 

.7019 

.7043* 

.  r. 

0. 

0. 

.  8  1  83 

-.0380 

.7798 

.7810* 

.  8 

0. 

0. 

.  9003 

- .  03  50 

.8097 

.  8070* 

0. 

12. 

12. 

.  7981 

-  .  0378 

.7004 

.7000 

.  K 

!  2. 

12. 

.  844  4 

-  .  0389 

.  8000 

.7942 

12 


fl 

II 

y! 

h 

0(  1  ) 

D(n*  )  »>( 

I  )M) (n~  ) 

Kxacl 

•10 

1.2397  0. 

0. 

0. 

.7758 

-.0272 

.7486 

.75* 

.  2 

0. 

0. 

.7804 

-.074 

.7530 

.7544* 

.  5 

0. 

0. 

.8093 

-.0285 

.7808 

.7818* 

.8 

0. 

o. 

.8968 

-.0272 

.8696 

.8683* 

0. 

20. 

20. 

.7893 

-.0278 

.7615 

.7564 

.8 

20. 

20. 

.8  352 

-.0290 

.8062 

.8038 

Next,  we  .study  the  accuracy  of  the  asymptotic  expression 
for  the  case  of  the  distribution  of  the  sample  correlation 
coefficient  rJ2=  r.  Wheny^y^  0.  The  distribution  of  r 
was  first  found  by  Fisher  (1915),  and  Hotelling  (1953)  has 
expressed  the  distribution  in  terms  of  the  hypergeometric 
function.  When  i i^O,  the  distribution  of  r  is  complicated 
and  the  cumulative  distribution  of  r  has  been  tabulated  by 
David  (1938).  Pillai  (HMD)  suggested  t  h*>  transformations 
g(r)  -  (r-p) /( 1-r c)  that  renders  the  distribution 

close  to  the  Student’s  t.  distribution.  Tn  Table  4  we 
compare  the  asymptotic  expressions  for 

Pri  r  £  ii|c  >  Y2)  and  Prfg(r)  u  I  p  .  y  ^  ,  y^  ) 

with  the  corresponding  values.  The  exact  values  with  *  are 
taken  from  the  tables  of  David  (1938)  and  the  remaining  exact 
values  are  obtained  by  simulation  with  5000  trials. 

TABLE  4 

Comparison  of  the  Asymptotic  Expression  with  Exact 
Expression  for  the  Distributions  of  Functions  of  the 
Sample  Correlation  Coefficient 

n  p  u  Yj  Y2  Stat.  Otl)  D(  n*  )  0(l)+0(n'~)  Kxacl 

40  .5  0.  0.  r  .5000  -.0142  .4858  .4856* 

l  7 


0. 

0. 

0. 

K(r> 

.5 

24 . 5 

12.25 

r 

.0 

24.5 

12.25 

g(r) 

.  5 

24.5 

24 . 5 

r 

0. 

24.5 

24.5 

f?(r  ) 

.35 

0. 

0. 

r 

-.  1818 

0. 

0. 

r(  r  ) 

.35 

24.5 

12.25 

r 

-.1818 

24.5 

12.25 

.35 

24 . 5 

24.5 

r 

-. iHI 8 

24.5 

24.5 

R(  r) 

.7 

0. 

0. 

r 

-.2273 

0. 

0. 

K(  r) 

.7 

12. 

6. 

r 

-.2273 

12. 

6. 

K(r) 

.7 

12. 

12. 

r 

- . 2273 

12. 

12. 

g(  r) 

.55 

0. 

0. 

r 

-.4464 

0. 

0. 

g(r) 

.55 

12. 

6. 

r 

-.4464 

12. 

6. 

g(  r) 

.55 

12. 

12. 

r 

■ifiiiTia 

1 

.4856 

.8748 

.0106 

.8854 

.8818 

.8948 

-.0088 

.8860 

.9235 

.0123 

.9358 

.9370 

.9438 

- . 0077 

,9361 

.0808 

.0156 

.0964 

.0966 

.  1016 

-.0063 

.0952 

.4486 

- . 0064 

.4422 

.4486 

.4491 

- . 0068 

.4422 

.  5568 

-.0044 

.  5524 

.  5586 

.  5573 

-.0049 

.5524 

.  0868 

.  0349 

.  1217 

.  1183 

.1328 

-.0175 

.1152 

.8341 

.0208 

.8549 

.8434 

.8800 

-.022 

.8576 

.9105 

.0295 

.9400 

.9370 

.  9600 

-.0180 

.9420 

.0003 

.0022 

.  0025 

.0098 

.0144 

- . 0030 

.0114 

.3870 

- . 0069 

.3802 

.3680 

.3924 

-.0124 

.3800 

.5534  -.0057  .5478 


4464  12 


12 


B(r )  . 5547 


0069  .5477 


5388 


i  X 


4.  Asymptotic  Distributions  of  Functions  of  the  Elements 
of  the  Sample  Covariance  Matrix  lor  Nonnormal 

Populations 

Let  X  =  t , . . . , X  3  where  the  p*l  random  vectors 

1 X  are  distributed  independently.  Also,  let 

S  =  XX'  =  (S^j)  be  *bo  samPle  sums  of  squares  and  cross- 

products  matrix  such  that  E(S/n)  =  £2  =  foi,  . ) ,  s .  .  =  S  /n  and 

i.l  ’  ij  ij' 

Y  -  •*(  I  -  a)‘  <yu> 

In  addition,  let  the  functions 


T^S/n)  Ti(sir...  .spp,s12,  .  .  .  *  sip  >  s23  >  •  •  •  *sp-ifP> 


be  analytic  in  the  neighborhood  of  fi  for  i=l,2,...,k.  The 
Taylor  expansion  of  Lj  =  Sn  (T^S/n)  -  T±( S2) )  about  ft  is 

Li  35  I  aj1j  yj  1  +  p  l  £  a(l^  y  y 

1  JA»J2  JlJ2  3l°2  2 ^  jl’J2  <1 3 * .14  J*J2  Js'14  J1J2  •13J< 


+  0(n_1) . 


Hence 


Uj  «=  E(L 


E  l  ai1i  -i  i  + 

1  2/n  jrj2  j3><j4  3iJ2  hU  1  *  6  4 


“ij  " 


E(L.L  )  =  l  l  »(i)  -<J> 

.1  J  J  i  . 


0(n  *) 


-1, 


^  =  jxfja  ^  “*1*2  *^4  +  0(0  } 


Uijt“  E(LiLjLt)s 


l  l  l  a  1 1  j  a ^ ^ j  ic(j1j2, j3j4,j5 

1 1 ,  j2  J3  >  J4  j5»J6  J1J2  J3>14  J5J6  1  *  J  *  5 


+  0  (n"1) 


(4.1) 
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where 


<(  j-i  jo.jojL)  =  E(y  .  .  y.  .  ) 
y  &  o  **  J i J2  J3J4 


<( j, Jo, J0J4. Jej«)  =  E(y .  .  y  .  y.  .  ). 

1  *  J  4  D  b  J1J2  J3J4  J5°6 


Now,  let  k .  .  ....  denote  the  cumulant  of  order  u  for 
J1J2  Ju 


Xi  where  . ju  =  1 . p.  Then 


Kjij2‘i3^4  +  <^1K^2'^3^4  +  (KjlJ3Kj2J: 


(4.2) 


J1J4  ^2^3^  ^  ^1^3  323‘ 


jjj2i ^3^4 , 


■  i  1C .  4-  £  K . K  .  +  £  (C  .  .  .  K  ,  .  . 

/n  J1J2<:!3J4J5J6  3  1323334  J5  J6  J1J2J3  J4'15J< 


J.  £. _ _  J.  A _ *  v _ _ _ 

+£  K  .  .  .  K .  +  £  K .  .  .  .  K .  K .  +  J  K .  .  .  K .  .  K . 

J1J2,13J5  J4°6  ■11J2J3J5  J4  J6  '  J1J2J3  J4J5  Jl 


(4.3) 


O _ _  _  O  _ _ .. 1 1  _ _ 

+ 1  K  .  .  .  K  .  K  .  ~7“  +  }'  <7  K  .  .  K  .  .  +  l  K  ,  ic  .  .  K  .  K  . 
J1J3J5  3 2  J4  ll6  J1J4  J2J5  J3J6  J1J4  J2J5  J3 


The  expressions  under  " - 

over  n  samples.  For  instance: 


represent  the  average  values 


n  k  .  T~7—  =  l  K(i> 

J1J3  J2J4  i=l  J1J3  J2J4 


The  summations  in  Eqs.  (4.2)  and  (4.3)  are  over  the  possible 
ways  of  grouping  the  subscripts  and  the  number  of  terms  re¬ 
sulting  is  written  over  jt.  fiqs .  (4.2)  and  (4.3)  coincide  with 
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the  expressions  of  Kaplan  (1952)  when  Xlf...X  are  identically 
distributed . 

So  the  cumulants  of  I,  =  (L  ,  .  .  .  ,L,  )  are 

~  jL  k 


*i =  “i 


‘ij =  “ij +  °<"  > 


(4.4) 


*ij£’“ijA  -  (V\j£+Uj“i£  +  U£Uij)+0(n”  > 


for  i , j , £=1 , . . . , k .  An  equivalent  equation  (see  Kendall 
and  Stuart  (1961)) 


=  ~  K(j1J2,  j3j4)tc(  j5J6,  j7jg)  +  °(n_1)  (4.5) 


is  used  in  calculating  «^u^,  n.u.f,  «}u.  .  in  Eq.  (4.4). 
The  approximated  characteristic  function  of  L  is 


E[exp(  i  t '  L)  ]  =  exp(  -  ^  t  Q  t  +  i  }'  t .  K. 

~  ~  *  '  1=1  3  1 


+  i3  y  U— r^r  t.t.  .t..K.  ..)  +  0(n'1) 

''  -6(1,  3  .  *)  1  3  l  ijl' 

*  i 


1 

exp(-  i  t  Q  t)  (1  +  i  l  t .  K, 
*  ~  i=]  11 


(4.6) 


*  ^  ,1  WLTJn  ‘i'j'.'y.1  +  0<"'1) 
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n- 


5.  Appl icat  ions  of  the;  Distributions  of  Functions  of  the 
Elements  of  the  Covariance  Matrix  and 
Correlation  Matrix 

In  this  section,  wo  discuss  some  applications  of  the 
results  of  Section  2  in  simultaneous  tests  of  hypotheses 
o.n  the  elements  of  the  covariance  matrix  and  correlation 
matrix. 

M-°'  "i.t:  "ij-’ou-  aU:  ‘Wj  •  "*J:  °ij'POij 

and  A,  p.  .4  i>  We  will  first  discuss  the  problem  of 
i.1  U  Oij 

testing  the  hypotheses  H. ^  simultaneously  against  the 
alternative  hypotheses  A^.  In  this  case,  the  hypothesis 
Hjj  is  accepted  if 

h.  <  s,  c  .5  al  *  for  i  <  j 


l  *  ‘  IJ  Oi.1 


*»ii 

£  0--  1  «2 

Oil 


for  i , j=l ,2 , . . . ,p ,  where 


(5.1) 


‘ii 


. ”• 

i“l . .  I  H  j  1  =  (  1  -  a) 

and  H.  =  D  II.  .  *  For  practical  purposes,  we  may  choose  the 
1  i< j  ' J 

constants  aj  and  b^,  si 2  an(*  ^2  suc^  that  a^=  -bj  and 

a^*  1/bg.  We  can  propose  similar  procedures  for  testing  the 

hypotheses  against  one-sided  alternatives. 
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(5.2) 


Next,  consider  the  problem  of  testing  the  hypotheses 

H.^ . Hpp  simultaneously  against  . A  .  In 

case,  we  accept  if 


s .  . 

i  <  11  <  n 

}2  —  a  —  a2 
*  Oil  * 


and  reject  it  otherwise  where 

Sj 


Ptb2  1 


ii 


-  <  a9;  i=l . p|  oH..]  =  (1-a).  (5.3 

Oii  “  *  i=l  11 


Also,  consider  the  problem  of  testing  the  hypotheses 

H*j(i  j=l,2,...,p)  simultaneously  against  the  alternatives 

A*  .  where  Ilf  .  :  p .  .  =  prt .  .  and  A*  .  :  p.  .  4  n_  •  •  •  In  this 

ij  ij  ij  Oij  i,i  1  xj  T  ^Oij 

case,  we  accept  or  reject  according  as 


C 

ct 


where 

Prr(rij-p0i j )2  <  ca*  i <j  =  l , 2 . p!H*J  =  (1-ot) 

and  H*  =  n  Ilf  . . 
i  <  j 

The  results  of  Section  2  are  useful  in  computing 
approximate  values  of  the  critical  values  associated  with 
the  above  tests.  The  results  of  Section  2  are  also  useful 
in  finding  approximate  critical  values  of  the  tests  of 
Krishnaiah  (1975)  for  testing  the  hypothesis  an'*,*‘”0pp 


2 


r 


against  different  alternatives  when  the  correlation  matrix 

is  known  as  well  as  the  procedure  of  Roy  and  Bargmann  (1958) 

for  testing  the  hypothesis  that  the  covariance  matrix  of 

multivariate  normal  is  diagonal. 

In  the  appl icat ions  discussed  above,  we  assumed  that  the 

matrix  S  *  (S.  .)  is  the  central  Wishart  matrix.  But 

i.J 

situations  arise  where  i he  model  itself  is  not  correct.  For 
example,  if  we  assume  that  X|,.--,\  are  distributed  inde¬ 
pendently  and  identically  as  multivariate  normal  with  a 

common  mean  vector  p  and  covariance  matrix  T. ,  then  R  is  the 

0  _  . 

central  Wishart  matrix  when  S’  =  )'  (X  .. -X.  )  (X  • -X  .  )  and 

n  j  =  l  ~J  ~J 

nX.  =  £  X..  But,  if  the  mean  vectors  are  given  by 

j  =  l 

F(N.)  =  p.  then  S  is  the  noncontra  1  Wishart  matrix.  So, 

-J  -J 

t  ho  results  given  in  Sect  ’on  2  for  the  case  of  ’the  noncentral 
Wishart  matrix  and  noncentral  correlation  matrix  are  useful 
in  studying  the  robustness  of  several  test  procedures  on  the 
elements  of  X  if  the  assumption  of  common  mean  vector  for 

Xll . Xn  is  viw,atod* 

Next,  consider  the  model. 

X-  .  =  u  .  +  6  . 

1.1  J  J 

X_ .  =  v  .  +  c  . 

2.1  J  J 

when  u.  =  a  +0v.,  (j=l,2 . n)  and  «,$  are  unknown  con- 

J  J 

stants.  Also,  we  assume  that  v.’s  are  distributed  inde- 

J 

pendently  and  identically  as  normal  with  a  common  mean  p 


and  variance  a  .  In  addition,  e . ' s  and  6.'s  are  distri- 

J 
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buted  as  normal  and 

E(c  )  =  E(S  )  =  0 

J  «J 

cov( e  . , 6  .  )  =  0,  cov(v.,e.)  =  cov(v.,6.)  =  0 
J  J  J  J  J  J 

Var(cj)  =  o2,  Var(6j)  =  a2. 


We  also  assume  that  the  random  vectors  (X  .  ,X,u.,v.,e,,<5.) 

"2 3  3  3  J  J 

are  distributed  independent  of  each  other  for  different 

2  2 

values  of  j.  When  X  =  /  ci£  is  known,  the  maximum  likeli¬ 

hood  estimate  of  0  is  known  (see  Kendall  and  Stuart  (1973), 
Chapter  29)  to  be 

a  _  (sn  S22)  +  U  Sj  1~X  S22)  +4AS]2^‘i 

3 - -t*— - 

^  12 


where  S  =  (  S.  .  )  =  ^  (  X.  -  X.  )  (  X.  -  X  )  '  ,  X =  ( X,  .  ,  X„  .  )  ,  and 

n  i  J  -i  =  1  -1  -J  ' 1 J  <2.1 

—  «J  1 

nX.  =  l  X  The  results  in  Section  2  of  this  paper  are  useful 
J 

in  obtaining  an  asymptotic  expression  for  the  distribution  of  0. 

Now,  let  y’  =  (Y.  Y,,),  ( j  =  1 , 2  ,  .  .  .  ,  n )  be  distributed 

~  J  -l  J  <2  J 

independently  as  a  bivariate  normal  with  mean  vector  (u1.,u0.) 

A  j  &  J 

2 

and  covariance  matrix  l.  -  o  T  where  I  i  s  an  identity  matrix. 

Also,  let  u  •  =  ex  +  0u  Then,  the  maximum  likelihood  esti- 

1 J  2  3 

mate  0  of  g  is  known  (e.g.,  see  Anderson  (1976))  to  be 


0  = 


S11-S22+«S11-S22>2+4S?2>J 


2  S 


12 


n  _  _  _  n 

where  S  =  (  S  )  =  l  (Y  .-Y.  )(Y  -Y.  )  ’  and  nY.=  [  Y, 

j_|  -J  ~J  ~  ~  j=l'J 
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Since  S  is  distributed  as  the  noncentral  Wishart  matrix, 

A 

an  asymptotic  expression  for  the  distribution  of  B  can  be 
detained  as  a  special  case  of  the  results  of  Section  2. 
Here,  we  note  that  Kunimoto  (1980)  has  recently  obtained  an 

A 

asymptotic  expression  for  the  distribution  of 
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